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ABSTRACT 

We have monitored the BL Lac object S5 0716+714 in five intermediate op- 
tical wavebands from 2004 September to 2011 April. Here we present the data 
that include 8661 measurements. It represents one of the largest databases ob- 
tained for an object at optical domain. A simple analysis of the data indicates 
that the object was active in most time, and intraday variability was frequently 
observed. In total, the object varied by 2.614 magnitudes in the i band. Strong 
bluer-when-brighter chromatism was observed on long, intermediate, and short 
timescales. 

Subject headings: BL lacertae Object: individual (S5 0716+714) - galaxies: active - 
galaxies: photometry 
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Introduction 



Blazars constitute the most variable subclass of active galactic nuclei (AGNs). 
Depending on whether or not showing strong emission lines in spectra. Blazar is 
divided into flat-spectrum radio quasars (FSRQs) and BL Lacertae (BL Lac) objects. 
BL Lac objects are characterized by non-thermal continuum emission across the whole 



electr omagnetic spectrum with absent or weak emission and absorption lines (IStickel et al. 



1993 1), variable and hi gh polarization ( jAngel fc Stockman 



1980 



Gabuzda et al. 



to gamma rays (IRavasio et al. 



Impey fc Neugebaur 



1988 



1989) , large amplitude and rapid variability at all wavelengths from radio 



radio components ( IDenn et al. 



2002 



2000). 



Bottcher et al. 



20031 ). and super luminal motion of 



S5 0716+714 is a distant BL Lac. In 1979, it was di scovered in a survey for sources 



19811 ). Because of the featureless 



with a 5 GHz flux greater than 1 Jy (jKuhr et al. 
spectrum and its strong optical polarization, it was identified as a BL Lacertae object by 



Biermann et al. 



(119811 ) . The redshift of S5 0716+714 was uncertain until Wagner et al. 
(1996) estimated a value bigger than 0.3. Afterwards, Nilsson et al. (2008) acquired a deep 
i-band image of this object and derived a redshift of 0.31 ± 0.08 by using the host galaxy as 
a standard candle. Most recently, Danforth et al. (2012) set an upper bound of z < 0.304 
with a confidence level of 90% for this object. This source is one of the most studied BL 
Lac objects, because it has high brightness and strong variability. The optical duty cycle of 
S5 0716 +714 is nearly unity, ind icating that the source is always in an active state in the 



visible (IWagner fe Witzel 



19951 ). Strong bluer-when- 



both internight and intranight variatio ns ( jWu et al. 



brigh t er co r relations were found for 



2005 



2007 



2012 



Poon et al. 



2009 



Hao et al. Il2010l ; I Chandra et al. 



201lh . 



This source has been intensively monitored by a number of authors. During a 4-week 
period of continuous monitoring, the source displayed in both optical and radio regimes a 
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transition between states of fast and slow variability with a change of the typical variability 



timescale from about 1 to about 7 days (jOuirrenbach et al 



19911). Wagner et al. (1996) 



investigated the rapid variations of this object in the radio, optical, ultraviolet, and X-ray 
regimes and found that it always keeps high amplitude change on the timescale of a few 

^ac to be ~0.4 mag 



Raiteri et al. 



((2003!) 



days. Sagar et al. (1999) showed an average V — R color of this BL 
in their one month long BVRI optical monitoring campaign in 1994. 
reported that the long-term optical brightness variations of this source appear to have a 
characteristic timescale of 3.3yr and four major optical outbursts were observed at the 
beginning of 1995, in late 1997, at the end of 2000, and in fall 2001. In particular, an 
exceptional brightening of 2.3 mag in 9 days was detected in the R band on 2000 October 
30. Color analysis on the optical light curves reveals only a weak general correlation 
between the color index and the source brightness. Recently, Poon et al. (2009) monitored 
the BL Lac object S5 0716+714 in the optical band during 2008 October and December 
and 2009 February with a best temporal resolution of about 5 minutes in the BVRI bands. 
Typical timescales of microvariability range from 2 to 8 hr. The overall V — R color index 
ranges from 0.37 to 0.59. Strong bluer-when-brighter chromatism was found on internight 
timescales. The overall variability amplitude decreases with decreasing frequency. 

We have monitored S5 0716+714 since 2004. Here we present the data during the 
period from 2004 to 2011. A simple analysis is performed and the results are described. 

This paper is organized as follows. The Observation and data analysis is described 
in Section 2. Section 3 presents the light curves. Section 4 shows the comparison result 
of i-data from us and i?-data from other authors. The relation of color and magnitude is 
described in Section 5. The conclusions are given in Section 6. 
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2. Observations and data analysis 

Our optical monitoring program of S5 0716+714 was carried out with the 60/90 
cm Schmidt telescope located at the Xinglong Station of the National Astronomical 
Observatories of China (NAOC). Prior to 2006, a Ford Aerospace 2048x2048 CCD camera 
was mounted at its main focus. The CCD has a pixel size of 15 fim, and its field of view 
is 58'x58', resulting in a resolution of l'!7 pixel -1 . At the beginning of 2006, the 2k CCD 
was replaced by a new 4096x4096 CCD. The field of view is now 96'x96', resulting in 
a resolution of 1'.'3 pixel -1 . The telescope is equipped with 15 intermediate-band filters, 
covering a wavelength range from 3000 to 10000 A. 

This paper includes data from 2004 September 10 to 2011 April 24. Excluding the 
nights with bad weather and those devoted to other targets, the actual number of nights 
for S5 0716+714 observations is 332. We used filters in e, i, and m bands to observe 
in 2004-2006, and then changed to the c, i, and o bands from 2006 December. The 
central wavelengths of the c, e, i, m, and o bands were 4210, 4920, 6660, 8020, and 9190 
A, respectively. The central wavelength of i band is similar to the R band. With the 
observational results of stars, t 



R = i + 0.1 flZhou et al. 



:he magnitudes in these two bands can be transformed with 
20031 ). Depending on the weather and seeing conditions, the 
exposure time of different bands range from 30s to 480s, and the exposures per night of 
different bands varies between 2 to 43. 

The data reduction procedure includes bias subtraction, flat-fielding, extraction of 
instrumental aperture magnitude, and flux calibration. We used differential photometry. 
For each frame, the instrumental magnitudes of the blazar and four comparison stars (See 
Fig.l) were extracted at first. The radii of the aperture and the sky annuli were adopted as 
3, 7, and 10 pixels, respectively. Then the brightness of the blazar was measured relative 
to the average brightness of the three reference stars 3, 4, and 5. Star 6 acted as a check 
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star, which has an apparently similar bright 



reference and check stars, see 



Howell et al. 



ness a s the blazar (for a reasonable selection of 



19881 ). The differential magnitude of star 6 is 



the difference between the magnitude of star 6 and the average magnitude of star 4 and 
5, so as to verify the stable fluxes of the four comparison stars, and to verify the accuracy 
of our measurements. The c, e, i, m, and o magnitudes of the 4 comparison stars were 
obtained by observing them and the standard star HD 19945 on a photometric night and 
are listed in Table 1. 



3. Light Curves 

The samples of observational log and results are given in Tables 2-6. The columns 
are observation date and time in universal time, Julian date, exposure time in second, 
magnitude and error of S5 0716+714, and differential magnitude of star 6 (its nightly 
averages were set to zero). Figure 2 shows the light curves of the overall monitoring period 
in the five bands. 

The source remained active during the whole monitoring period. The variation 
amplitudes of e, i, and m bands from 2004 September 10 to 2006 March 29 are 1.200, 
1.156, and 1.127 mags, respectively, and the variation amplitudes of c, i, and o bands from 
2006 December 6 to 2011 April 24 are 2.763, 2.614, and 2.522 mags, respectively. The 
amplitude of variation tends to decrease with decreasing frequency. The light curves keep 
fluctuating during the whole monitoring period. The curves get the extreme bright value 
on 2004 September 10 (JD 2,453,259), 2007 October 20 (JD 2,454,394), 2008 April 22 (JD 
2,454,579), and 2010 September 28 (JD 2,455,468), and get the extreme dark value on 
2005 January 29 (JD 2,453,400), 2007 December 16 (JD 2,454,451), and 2011 April 24 (JD 
2,455,676). A faintest optical state was recorded on 2007 December 16 (JD 2,454,451), . 
Nilsson et al. (2008) acquired a deep i-band image of this object at that state and derived 



-7- 



a redshift of 0.31 ± 0.08 by using the host galaxy as a 'standard candle'. 



The BL Lac object S5 0716+714 is one of the brightest BL Lac objects noted for its 
microvariability. In order to confirm whether or not the object was variable in one day, 
a quantitative assessment was carried out. For each of the 233 nights with observational 



duration longer than 2 hours, we used a chi-square inspection to check whet 



rer there 



is int raday variability (IDV) (IPenston fc Cannon 



1970 



Kesteven et al. 



1976 



de Diego 



2010J). The chi-square value was compared with the critical value at the 95% confidence 
level. If the former was greater than the latter, the null hypothesis that there was no 
variability was rejected. As a result, 138 nights (62% of 233) with IDV were identified. The 
observed fastest variation of S5 0716+714 varied by 0.117 mags in 1.1 hrs in the c band 
on 2011 March 5 (JD 2,455,626), as indicated by the bottom line with arrows in Figure 3. 
The similar work was made by Villata et al. (2000). They found that the source exhibited 
strong variability with similar trend but different amplitudes in all bands, and noted the 
monotonic brightness increase in B band for about 130 minutes. The steepest (linear) part 
has a rising rate of 0.002 mag per minute and a duration of about 45 minutes. 



4. Compare with Other Data 

In Section 2, a transforming formula between the i and R magnitudes was mentioned. 
This formula was derived from the observational result of stars (Zhou et al. 2003). However, 
the spectral shape of blazars is quite different from that of stars. So the transforming 
formula may be different for blazars. Therefore, we made a comparison between our z-band 
data and the i?-band data of other authors in order to find an empirical relation between 
them. Villata et al. (2000) and Poon et al. (2009) have made intensive monitoring on the 
same object. Their data were adopted and matched in time with ours with a threshold of 
less than 0.02 days (or 28.8 minutes). As the results, 24 R-i matches were found for Villata 
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et al. and our data, and 97 R-i matches were found for Poon et al. and our data. The 
average time differences are 0.0038 and 0.0009 days for the 20 and 97 matches, respectively. 
Two R-i diagrams were plotted in the top and middle panels of Figure 4. Two linear 
regressions give the transforming formulae as R = (0.964 ± 0.015) x i + (0.279 ± 0.187) 
and R = (0.976 ± 0.009) x i + (0.059 ± 0.123), respectively. If all matches are plotted 
together, as shown in the bottom panel of Figure 4, the linear regression gives the formula 
as R = (0.897 ± 0.004) x i + (1.127 ± 0.048). 



5. Relation of Color and Magnitude 



The long-term color behavior of S5 0716+714 was studied based on our data. For the 
2004-2006 data, the e — m color was calculated and plotted vs. the e magnitude in the top 
panel of Figure 5. For the data after 2006 December, the c — o color was calculated and 
plotted vs. the c magnitude in the central panel of Figure 5. The bottom panels illustrate 
how the c — o color and c magnitude changed with time. Despite the discontinuity in the 
top panel and significant scatter in the central panel, there is an overall bluer-when-brighter 
chromatism in both panels. 

In order to investigate the color behavior of S5 0716+714 on the intermediate timescale, 
three episodes in our monitoring are isolated. They are from JDs 2,454,101 to 2,454,115, 
from JDs 2,454,429 to 2,454,463, and from JDs 2,455,597 to 2,455,629. These three episodes 
lasted from two weeks to more than one month. In these periods, we have relatively 
continuous monitoring, and the object showed significant variations. The light curves and 
the corresponding color-magnitude diagrams are plotted in Figure 6 for the three episodes. 
There are strong color-magnitude correlations. The correlation coefficients are 0.89, 0.78, 
and 0.82, respectively. The bluer- when-brighter chromatism on in termediate timescale is 



found by other authors (e.g., 



Villata et al. 



2000 



Wu et al. 



20071 ). On intraday timescales, 
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the object also displays strong bluer-when-brighter chromatism. Some examples are given 
in Figure 7. 

Depending on the balance between escape, acceleration, and cooling of the electrons 
with different ener gy, either soft (low energy) or hard (high energy) lags are expected 



( jKirk et al. 



1998). This will lead to a loop-like path of the blazars state in a color- 



magnitude (or spectral index-flux) diagram. The direction of this spectral hysteresis can 
be either clockwise or anticlockwise. It depends on the relative position between the 
observing frequency and the peak frequency of the synchrotron component in the SED of 
the blazar as well as the relative values of th e acceleration, cooling and escape timescales 



( iChiaberge fc Ghisellini 



1999: 



Dermer 



19991 ). 



In our monitoring, we got an inconspicuous loop flare in internight timescale and an 
obvious loop flare in intranight timescale, which are shown in Figure 8. The upper-left 
panel is the light curve of JDs 2,454,390 ~ 2,454,398, which shows a flare in do bands. 
The corresponding color-magnitude diagram is displayed in the upper-right panel, in which 
the numbers denote time sequence. There were not significant variations during these day, 
so we averaged the data in day to decrease system error. The points spread as diagonal 
distribution in color-magnitude diagram. The loop flare is far less obvious. The result of 
intranight timescale, JD 2,455,621, is shown in the bottom of Figure 7. The numbers denote 
the time sequence. An anticlockwise loop can be seen in the lower-right panel. As Kirk et al. 
(1998) supposed, if the loop is traced anticlockwise, there might be a flare propagating from 
lower to higher energy, as particles are gradually accelerated into the radiating window. 
The frequency of c and o band are 7.13 x 10 14 Hz and 3.26 x 10 14 Hz. In 1999, Giommi et 
al. found that the frequency is between 10 14 and 10 15 Hz. Our anticlockwise loop imply 
that the peak frequency of short band in the SED of the blazar should be higher than 
7.13 x 10 14 Hz. 
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6. Conclusions 

We have monitored the BL Lac object S5 0716+714 in five intermediate optical 
wavebands from 2004 September to 2011 April by the 60/90 cm Schmidt telescope located 
at the Xinglong Station of the National Astronomical Observatories of China (NAOC). We 
collected 8661 data points with error less than 0.05 mags. It represents one of the largest 
databases obtained for an object at optical domain and can be used to study both the long- 
and short-term flux and spectral variability of this object. It can also be correlated with the 
data in the radio, X-ray, or gamma-ray wavelengths in order to investigate the broad-band 
behavior of this object. A simple analysis of the data indicates that the object was active 
in most time. The overall amplitudes of e, i, and m bands from 2004 September 10 to 2006 
March 29 are 1.200, 1.156, and 1.127 mags, respectively, and the overall amplitudes of c, 
i, and o bands from 2006 December 6 to 2011 April 24 are 2.763, 2.592, and 2.522 mags, 
respectively. The amplitude of variation tends to decrease with decreasing frequency. The 
bluer-when-brighter phenomenon is effectively confirmed on long, intermediate, and short 
timescales. It is an important support to the shock-in-jet model in which shocks propagate 
down the relativistic jet, accelerating particle s and /or compressing magnetic fields, leadi ng 



to the observed flux and spectral variability (IMarscher fc Gear 



1985 



Qian et al. 



19911). 



There were 138 nights of IDV captured during the whole monitoring period. 
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7:23:00 7:22:30 7:22:00 7:21:30 7:21:00 



Right ascension 



Fig. 1. — Finding chart of S5 0716+714 and the 4 comparison stars taken with the 60/90 
Schmidt telescope and filter i on 2011 April 12 (JD 2,455,664). The size is 12' x 12' (or 
512x512 in pixels). 
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Fig. 2. — Light curves of S5 0716+714 in the c, e, i, m, and o bands. 
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Fig. 3. — The observed fastest variation of S5 0716+714. The object varied by 0.117 mags 
in 1.1 hrs in the c band, as indicated by the bottom line with arrows. 
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Fig. 5. — Color-magnitude diagrams of em bands (top) and co bands (middle) and the color 
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Fig. 6. — Light curves (left) and Color-magnitude diagrams (right) of intermediate 
timescales. The line in the right figures is the linear fit to the data. 
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Fig. 7. — Light curves (left) and Color- magnitude diagrams (right) of IDV data. The line in 
the right figures is the linear fit to the data. 
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Fig. 8. — Light curves (left) and color- magnitude diagrams (right) for the internight (upper) 
and intranight (lower) timescales. The numbers in the right panels denote the time sequence. 
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Table 1. Magnitudes of 4 Comparison Stars. 



Star 


c 


e 


i 


m 





3 


13.033 


12.706 


12.271 


12.190 


12.178 


4 


13.574 


13.412 


13.098 


13.050 


13.060 


5 


14.126 


13.794 


13.363 


13.264 


13.223 


6 


14.195 


13.890 


13.442 


13.370 


13.329 



Table 2. Data of c band 



Date(UT) 


Time 


Julian Date 


Exp 


c 


c err 


dfmag 


2006 12 06 


18:43:43.0 


2454076.28036 


300 


14.649 


0.016 


0.002 


2006 12 06 


18:58:33.0 


2454076.29066 


300 


14.636 


0.020 


-0.018 


2006 12 06 


19:12:43.0 


2454076.30050 


300 


14.639 


0.014 


0.011 


2006 12 06 


19:27:06.0 


2454076.31049 


300 


14.660 


0.014 


0.004 


2006 12 06 


19:41:42.0 


2454076.32063 


300 


14.612 


0.014 


-0.003 



Note. - Table [2] is published in its entirety in the electronic edition of 
the The Astrophysical Journal Supplement. A portion is shown here for 
guidance regarding its form and content. 



-26 - 



Table 3. Data of e band 



Date(UT) 


Time 


Julian Date 


Exp 


e 


c err 


dfmag 


2005 01 28 


15:57:45.0 


2453399.16510 


300 


14.216 


0.027 


-0.009 


2005 01 28 


16:07:44.0 


2453399.17204 


300 


14.212 


0.026 


-0.028 


2005 01 28 


16:17:43.0 


2453399.17897 


300 


14.186 


0.023 


0.015 


2005 01 28 


16:27:57.0 


2453399.18608 


300 


14.239 


0.029 


-0.015 


2005 01 28 


16:37:53.0 


2453399.19297 


300 


14.242 


0.026 


0.011 



Note. - Table [3] is published in its entirety in the electronic edition of 
the The Astrophysical Journal Supplement. A portion is shown here for 
guidance regarding its form and content. 



-27- 



Table 4. Data of i band 



Date(UT) 


Time 


Julian Date 


Exp 


% 


t err 


dfmag 


2004 12 20 


17:45:00.0 


2453360.23958 


240 


12.739 


0.011 


-0.012 


2004 12 20 


17:51:11.0 


2453360.24388 


240 


12.745 


0.010 


0.010 


2004 12 20 


17:55:40.0 


2453360.24699 


240 


12.717 


0.010 


-0.026 


2004 12 20 


18:00:06.0 


2453360.25007 


240 


12.748 


0.011 


0.032 


2004 12 20 


18:04:36.0 


2453360.25319 


240 


12.733 


0.014 


-0.008 



Note. - Table H] is published in its entirety in the electronic edition of 
the The Astrophysical Journal Supplement. A portion is shown here for 
guidance regarding its form and content. 



-28 - 



Table 5. Data of m band 



Date(UT) 


Time 


Julian Date 


Exp 


m 


Tfl err 


dfmag 


2005 12 21 


18:46:03.0 


2453726.28198 


300 


12.549 


0.007 


-0.019 


2005 12 21 


19:02:20.0 


2453726.29329 


300 


12.550 


0.007 


-0.013 


2005 12 21 


19:18:14.0 


2453726.30433 


300 


12.551 


0.007 


0.012 


2005 12 21 


19:34:46.0 


2453726.31581 


300 


12.545 


0.008 


0.007 


2005 12 21 


19:50:41.0 


2453726.32686 


300 


12.537 


0.008 


-0.002 



Note. - Table [5] is published in its entirety in the electronic edition of 
the The Astrophysical Journal Supplement. A portion is shown here for 
guidance regarding its form and content. 



-29 - 



Table 6. Data of o band 



Date(UT) 


Time 


Julian Date 


Exp 





o err 


dfmag 


2006 12 06 


18:51:49.0 


2454076.28598 


300 


13.580 


0.018 


0.004 


2006 12 06 


19:06:42.0 


2454076.29632 


300 


13.555 


0.017 


0.007 


2006 12 06 


19:21:04.0 


2454076.30630 


300 


13.557 


0.016 


0.008 


2006 12 06 


19:35:38.0 


2454076.31641 


300 


13.562 


0.016 


-0.005 


2006 12 06 


19:50:16.0 


2454076.32657 


300 


13.532 


0.016 


0.004 



Note. - Table [6] is published in its entirety in the electronic edition of 
the The Astrophysical Journal Supplement. A portion is shown here for 
guidance regarding its form and content. 



